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Abstract: We present numerical experiments on dynamical friction in pairs of elliptical galaxies of 
unequal mass. We confirm that the self-gravity of the response is not important and show the drastic 
effect of the deformability of the companion which reduces the decay time by more than a factor 
2. Almost the same amount of orbital energy is dissipated within the satellite than within the large 
galaxy. Finally, we discuss the importance of distant encounters for the dynamical evolution of systems 
of galaxies. 

1. Introduction 

This communication presents the first results from a N-body study of pairs of elliptical galaxies, they 
concern the investigation of the effect of various parameters on the dynamical friction in unequal mass 
pairs. 

The important problem of the orbital decay of a “small” satellite onto a massive galaxy has been 
the object of a very open controversy over the last decade. The debate include the question of the 
globality of the response of the large galaxy and that of the importance of its self-gravity (Lin and 
Tremaine, 1983; White, 1983; Bontekoe and van Albada, 1987 (—BvA); Zaritsky and White, 1988 
(=ZW); Hernquist and Weinberg, 1989). Although the most recent publications agree on the merging 
time and, to a very satisfying accuracy, on the shape of the decay curve (separation versus time), the 
actual interpretation of the physical phenomena is still discussed. According to BvA, the dynamical 
friction can be characterized as a local process and the global response is not dynamically important. 
ZW agree about this last point, but claim for the globality of the nature of the dynamical friction 
although it can be described by a local formula. This discussion obviously suffers from the lack of 
clear and well accepted definitions of the terms. In this respect, the recent claim by Weinberg (1989) 
of the strong importance of the self-gravity, in the sense it increases the decay time by a factor 2, is 
very confusing. Indeed, Weinberg includes in the self-gravity the motion of the center of gravity of 
the large galaxy while the other only consider the effect of the deformations on the internal potential. 

In the present work we concentrate on the case of the friction at distance (when the two galaxies are 
well separated) which seems important for describing the earliest stages of the decay since the loose 
random encounters are strongly more usual in a cluster environment than the close ones. In addition, 
unlike previous work except Borne (1984), our companion is represented by a system of particles, and 
thus has the capability to absorb orbital energy as the large galaxy does. 

The next part presents the experimental setup. The incidence of two physical ingredients is studied 
in 3 and 4: The effects of self-gravity and freedom degrees in the satellite. Then the importance of 
the friction at distance is investigated. 

2. Experimental setup 

We started from a self-gravitating system of two spherical galaxies since the a priori allowed 
approximations were not obvious. Since the internal potential of the galaxies is computed on a 
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cubic grid of finite resolution, we added a rigid core in order to simulate the high central density peak 
characteristic of elliptical galaxies without using a too large number of particles. The equations of 
motion are: 


Xcl = mi “ ®< 5 l )Pkc + m 2 ~ X ^)Ple + m c2(x c 2 ~ X c i )pj (la) 

*1 12 

Xji = mi - Zji)Pkj + m 2 ^(*.12 - xji)pjj 3 * * * * * * 10 + • • • 

fcl 12 

+ m c i(* c i - Xji )p-fi + m c2 (x c2 - xji)p~ 2 j (16) 

Where we adopted the notations from White (1983), the indices 1 and 2 refets to each galaxy and 
should be inverted to obtain the accelerations of the second galaxy, mi and m 2 are the masses of the 
particles and m c \ and m C 2 those of the cores. 

For technical reasons the interaction between the particles of the two galaxies (second term of lb) 
is approximated by a multipolar expansion (White, 1983). This approximation is a crucial point, 
but since for the experiments presented here we have made further simplifications which avoid the 
difficulty, we will not address the point. 

Two cases are considered: (1) a rigid satellite. (2) a satellite represented by N-particles, without 
self-gravity. 

In case (1) the second terms of Eq. la and b are ignored, and all the mass of the satellite is concentrated 
in the rigid potential. This makes the code formally equivalent to the classical problem of the decay 
of a satellite. The energy, angular and linear momentum are conserved. 

In case (2) the two first terms of Eq. lb are replaced by the interaction with fixed potential: 

2 m i (xji) + ^m 2 V $2 (xji) 

kl 12 

The problem is then that considered by Borne (1984), and it does not formally conserve energy. 

In our code, running on the CCVR CRAY 2, each galaxy is made by 20000 particles, and the self- 
gravitating potential is computed on a cubic 60x60x60 points grid. Forces are derived and interpolated 
from this grid. In the case (2) only the symmetric spherical part of the initial (relaxed) potential is 

taken into account. 

The initial model has a central rigid Plummer containing 1/3 of the total mass and the particles axe 
distributed according to a truncated 1/r 2 density profile. The velocities follow an isotropic maxwellian 
distribution truncated at the escape velocity. Then, each galaxy is virialized separately, and the final 
density profile appears reasonably realistic. An important step, as stressed by Borne (1984), is the 
initial relaxation in the tidal field of the second galaxy before starting the interaction. We chose to not 
relax the initial galaxies in the actual initial tidal field, since such procedure is likely to exaggerate the 
phase coupling and the amplitude of the tidal bulges. Instead, we relax the galaxies in a spherically 

averaged tidal field. This has the advantage to avoid the energy shock responsible for strong transient 
effects (affecting the decay time by a factor two in Borne’s experiments). 

The mass ratio between the two galaxies lies in the range 3 to 10, typical of the distorted pairs of 
ellipticals presented by Prugniel and Davoust (this meeting). The mass of the primary galaxy is 

10 U M & and its radius 30 kpc, the radius of the companion is 10 kpc. 

In all the experiments presented here the companion was launched on a circular orbit of radius in the 
range 30 - 60 kpc. 
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As a check of our procedure we run an experiment comparable to the standard experiment defined by 
White (1983). The decay curves are compared in Fig. 1. The agreement on the merging time is very 
good, but there are some difference in the shape of the decay curves. Although we did not check it, 
we believe the differences are due to the different mass profiles between these experiments. 


Figure 1: Comparison with White 

(1983) experiment: decay of a point- 

mass satellite onto a large galaxy from 

a circular grazing orbit. The mass ratio ^ 

is 10 and the mass of the large galaxy M 

is 1O 12 M0. The broken line is White’s ^ 

data and the full line is our experiment. O 
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3. Effect of the self-gravity of the response 

The first point investigated is the role of the self-gravity of the response. This is done by comparing 
experiments with the version (1) of the code with other done with version (2) but with a rigid satellite. 
Fig. 2 presents the decay curves for these experiments. The differences are not larger than the 
dispersion between experiments started from various initial conditions. This confirms the results of 
BvA and ZW and shows the consistency of the two versions of the code. 


Figure 2: Effect of the self-gravity of 
the response of the large galaxy. Decay 
of a point mass satellite; same experi- 
mental setup as for Fig. 1. The full line ^ 
includes the self-gravity of the defer- 
mation (case 1) while the dashed line g 
does not (case 2). The difference is not .2 
significant. 
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4. Effect of the deformability of the companion 

The second point is the analysis of the effect of the ability of the satellite to absorb energy. Fig. 3 
shows the results. The effect is considerable: the merging time is reduced by a factor 2.5 and scanning 
the range of realistic values for the size of the satellite it is always decreased by a factor 2 to 3. 


Figure 3: Effect of the deformabil- 
ity of the satellite. Decay of an ex- 
tended satellite (radius 10 kpc) start- 
ing on a circular orbit. The masses are 
10 12 Mq and 10 11 Mq, the radius of the 
large galaxy is 30 kpc. The full liue is 
for a rigid satellite and in the experi- 
ment represented by a dashed line both 
galaxies are made with 20000 particles. 
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5. Importance of the friction at distance 


Clearly, the classical formula for the dynamical friction does not apply when the galaxies are well 
separated (since in this model the drag is proportional to the local density). However, it is important 
to know the relative importance of loose encounter over close encounter for what concern the evolution 
of groups of galaxies. The alternative approach is the impulsive approximation (Spitzer, 1958), but it 
is not suitable to our problem where the relative velocity is of the same order or even smaller than the 
typical internal velocities. However, Dekel et al. (1980) have shbwn that the impulsive approximation 
still gives results compatible with numerical experiments for slow hyperbolic encounters. 

We have conducted experiments starting from distant circular orbits. The difficulty of these 
experiments with low friction is the high sensitivity to the noise due to the limited number of particles. 
With 20000 particles, we were able to run experiments up to an initial orbital radius of 60 kpc (with 
a rigid satellite). The initial decay rates ( dR/dt ) are compared with the prediction of the impulsive 
approximation on Fig. 4. The dashed line corresponds to the classical impulsive approximation 
(Spitzer, 1958; dR/dt oc ft -3 ) while the full line is the restriction to the drag due to the particles with 
a velocity smaller than the encounter velocity ( dR/dt oc .ft -9 / 2 ). 
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Figure 4: Dependence of the initial 
decay rate ( dR/dt ) on the radius of ^ 
the. orbit. The masses are 10 12 Mq and ® 
1O 11 M 0 and the satellite is rigid. The o 
dashed line corresponds to the classi- \ 
cal formula for impulsive approxima- ^ 
tion (Spitzer, 1958), in principle not 5* 
suitable to this case, and the full line g> 
is the friction computed under the im- % 
pulsive approximation restricted to the h 

low velocity particles. S' 

<0 
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6 . Conclusions 

We have confirmed that the self-gravity of the response is not important and revealed the drastic effect 
of the ability of the small companion to absorb orbital energy. The amount of energy absorbed by the 
satellite is comparable to that dissipated within the large galaxy. This effect is associated to the very 
strong deformation of the satellite illustrated in Fig. 5. The distortion pattern of the simulated pairs 
looks very much like the one observed in real pairs (e.g. Davoust and Prugniel, 1988 or Prugniel and 
Davoust, this meeting). As for the global-local character of the friction phenomenon, we investigate 
this question in a subsequent publication (Prugniel and Combes, 1990). 


Figure 5; Tidal distortion of the galax- 
ies. This plot shows the density of par- 
ticles in a logarithmic scale at an ad- 
vanced stage of one of our experiments. 
The initial separation was 40 kpc and 
the mass ratio is 10. The field repre- 
sented is 120 kpc wide and the actual 
separation is 10 kpc. The merging will 
occur in 5 x 10 7 years. Note the ex- 
tended diffuse tail behind the satellite 
and the characteristic pattern of distor- 
tion symmetric with respect to the pair 
center. 
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